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Abstract 

We prove that the gravitational binding energy J7 of a self gravitating system described by a 
mass density distribution p(x) admits an upper bound B[p{x)] given by a simple function of an 
appropriate, non-additive Tsallis' power-law entropic functional Sq evaluated on the density p. 
The density distributions that saturate the entropic bound have the form of isotropic g-Gaussian 
distributions. These maximizer distributions correspond to the Plummer density profile, well 
known in astrophysics. A heuristic scaling argument is advanced suggesting that the entropic 
bound B[p{x.)] is unique, in the sense that it is unlikely that exhaustive entropic upper bounds not 
based on the alluded Sq entropic measure exit. The present findings provide a new link between 
the physics of self gravitating systems, on the one hand, and the statistical formalism associated 
with non-additive, power-law entropic measures, on the other hand. 
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I. INTRODUCTION 



The application of non-additive, power law entropic functionals (and the concomitant 
g-maxent distributions) to the study of self-gravitating systems has attracted considerable 
interest in recent years . In particular, the above mentioned entropic measures and the 
associated g-maxent distributions were found to be relevant for the study of astrophysical 
self-gravitatin g sy stems such as galaxies and galaxy clusters 12|-|l8| and of other astrophys- 



ical problems 



19 



24| . such as the statistical characterization of the distribution of stellar 



rotational velocities in star clusters, where g-maxent distributions were successfully applied 
by Carvalho, do Nascimento, Silva and De Medeiros 
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Compelling evidence for a connection between the physics of self-gravitating systems, on 
the one hand, and the statistical formalism associated with Tsallis non-additive power law 
entropy Sg 
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-27j. on the other one, was provided by the long term, A^-body numerical 
simulations performed by Taruya and Sakagami These authors studied the behaviour 
of a stellar system confined within an adiabatic wall that evolves from initial conditions 
not described by the Boltzmann-Gibbs distribution and found that, before entering the 
gravothermally unstable regime, the evolving self-gravitating system can be described re- 
markably well by a sequence of polytropic distributions with a time dependent polytropic 
index. These polytropic distributions have precisely the form of Tsallis' g-maxent distri- 
butions maximizing an Sg entropy under the constraints imposed b y t he total mass and 
total energy of the self-gravitating system, as was first pointed out in 28(|. Taruya and Sak- 
agami also reported numerical evidence indicating that the system exhibits the same kind 
of behavior even if the simulations are performed without an outer boundary. 

Exact results concerning the equation of state within Tsallis thermostatistical framework 
of a two-dimensional gas of gravitationally interacting particles were obtained by Sau Fa and 
Lenzi |9|. The Jean's instability of a "gravitational plasma" endowed with a g-gaussian veloc- 
ity distribution was analyzed by Lima, Silva, and Santos [8']. The gravothermal catastrophe 
of self-gravitating systems characterized by polytropic [that is, g-maxent) distributions was 
investigated in detail by Taruya and Sakagami (see js], O] and references therein). Chavanis 
and Sire [lO| showed that the criterion for nonlinear dynamical stability in the case of spher- 
ical stelar systems governed by the Vlasov-Poisson equations exhibits some resemblances 
with a criterion of thermodynamical stability. On such basis, Tsallis' functional Sg can be 
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regarded as a useful if-function that allows for establishing a continuous connection between 
stellar polytropes, on the one hand, and isothermal stellar systems on the other one. An 
increasing body of evidence has been reported in recent years indicating that the polytropic, 
q-raaxent distributions are useful to describe dark mater halos of galaxies and galaxy clusters 



12 



20|. 



Given a mass density distribution p(x) it is intuitively clear that its gravitationally bind- 
ing energy Q tends to decrease as the system becomes less localized. Entropic functional 
provide natural quantitative measures of the degree of "delocalization" of the density p. 
Consequently, one may expect that there is some relation between the binding energy Q and 
some appropriate entropic functional evaluated on the density p. The aim of the present 
work is to prove that such a relation does indeed exist. 

II. TSALLIS STATISTICS AND GRAVITATION 

Some of the most important links between the physics of self-gravitating systems and 
the statistical formalism based upon Tsallis' entropy come from the study of self-gravitating 



systems described by the Vlasov-Poisson equations 



In order to analyze some important 



features of galaxies and other self-gravitating systems it is useful to regard them as non 



coUisional systems constituted by particles (the stars) interacting gravitationally j29|. 
Within this approach, each star moves under the effect of the average gravitational potential 
$(x) due to the remaining stars in the galaxy. The state of the system is then given by a 
distribution function F(x, v,t) in position-velocity space, where x denotes a star's position 
vector, V stands for its velocity vector and t is the time. The quantity F(x, v, t)(ix(iv gives 
the number of stars having position and velocity vectors respectively within the volume 
elements (ix and d\ in position and velocity spaces. 

The dynamical equation governing the evolution of collisionless system is the Vlasov 
equation 

dF 

— +vV,F-V,$- V,F = 0, (1) 
coupled to the Poisson equation 

V^* = AnGp, (2) 
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with p(x, t) = m J F(x, v,t)(iv, or equivalently to, 



$(x,t) = -Gmj l^2^^dx'dV, (3) 

which relates the distribution function -F(x, v, t) with the Newtonian gravitational potential 
$. Here G denotes Newton's gravitational constant and m is the mass of each star (we 
assume that all stars have the same mass). 

The Vlasov equation can also be cast under the guise {dF/dt)orh. = 0. This means that 
the distribution function F adopts a constant value along the orbit of any star moving in the 
gravitational potential $. It is important to stress that the non-collisional treatment based 
upon the Vlasov-Poisson system yields only an approximate description of the dynamics of 



the system. This approximation is, however, a very useful one for various stellar systems 29 |. 
A complete description at the Newtonian level of the evolution of a self-gravitating A^-particle 
system is given by the complete set of N coupled Newtonian differential equations of motion 
corresponding to the particles. This approach to the dynamics of self-gravitating systems 
is usually investigated by recourse to numerical A^-body simulations (see, for instance, 
The total mass. 



M = m I F(x,v)(ixc/v (4) 

the total energy. 



„ m I ^ , , Gm? f F(x, v)F(x', v') ,,,,,, 
E = — v2F(x,v)dxrfv / ^ ^ \ ' — -dxdvdx.' dv', (5) 



2 J ' ' 2 J |x-x 

and the total angular momentum 



L = my XX vF(x, v) (ixfiv (6) 

are preserved under the time evolution determined by the coupled Vlasov ([1]) and Poisson 
([2]) equations. The Vlasov-Poisson dynamical system ([UlS]) also admits a family of conserved 
quantities given by functionals of the form 

K[F] = J g{F)d^dv. (7) 
In other words, time dependent solutions -F(x, v, t) of the Vlasov equation ([1]) satisfy 
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When one considers different stationary solutions of the Vlaosv-Poisson system as possible 
candidates for describing galaxies or other astrophysical A^-body systems, it is natural to 
ask which is the distribution function -F(x, v) that maximizes Boltzmann entropy under the 
constraints given by the total mass and total energy of the system. Unfortunately, the 
formal solution of this variational problem is given by a distribution function (the isothermal 
sphere) exhibiting the physically unacceptable features of having infinite mass and infinite 
energy 29|. On the other hand, the optimization of the non additive Tsallis entropic measure 
Sq under the same constraints yields the stellar polytropic sphere distributions which, for a 



certain range of va 
finite total energy 



ues of the entropic parameter q, are characterized by a finite mass and a 



The stellar polytropic sphere distributions in position-velocity space are given by 



F(x, v) = F(e) = A{% - e)"-3/2 g < $o (9) 
= e > $0, 

where 

e=^v2 + $(x), (10) 

denotes the total energy (per unit mass) of a single, individual star, and A, $o, and n 
(usually referred to as the polytropic index) are appropriate constants. After making the 
identification of parameters determined by 

^ n-l. (11) 



1-g 2 

it can by immediately seen that the stellar polytropic sphere distributions have the g-MaxEnt 
form. That is, they arise from the maximization of the Tsallis' functional 

Sq[F] = -1_ (^1 _ y F(x,v)^rfxdv^ , (12) 

under the constraints given by the total mass (jlj) and total energy ([5]). Furthermore, the 
velocity distribution corresponding to a given space location is a g-Gaussian. The isothermal 



sphere distribution is recovered in the hmit n — t- oo (corresponding to g — )■ 1). Another 
important special case is obtained for n = 5, which yields the Schuster sphere distribution 



29|. 

For n < 5 (and q < 7/9) the polytropic spheres are systems of finite total mass and energy. 
The cut-off exhibited by the polytropic distribution corresponds to what is known in the 
field of non extensive thermostatistics as "Tsallis cut-off prescription" . The physical origin 



of this feature of the stellar polytropic distributions is given 
velocity corresponding to each value of the radial coordinate r 



jy the gravitational escape 



29|. 



III. AN ENTROPIC UPPER-BOUND ON GRAVITATIONAL BINDING EN- 
ERGY 

The gravitational binding energy of a mass distribution characterized by the density 
p(x) is given by 

n = G [ 4^^dx dx'. (13) 
J |x X I 

The system's total gravitational potential energy is equal to —1]. In the case of a self- 
gravitating system constituted by particles of mass m and described by a space-velocity 
distribution -F(x, v) we have p(x) = J -F(x, v)(iv and 

n=^/£^^)^rfx<ivrfxVV. (14) 

Z J |X X I 

It is intuitively clear that the gravitational binding energy decreases as the system's 
mass distribution p(x) becomes more delocalized. Quantitative measures for this degree of 
delocalization are provided by entropic measures evaluated on a probability density defined 

as 

where M = J p(x)(ix is the total mass of the system. Note that /(x) is normalized to 
1. These considerations suggest that there may be a connection between the gravitational 
binding energy Q, on the one hand, and an appropriate entropic measure evaluated on the 
probability density ]5(x). Here we are going to prove that, indeed, such a connection exists 
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and can be established as an entropic upper bound on the binding energy Q. We are going 
to obtain an upper bound to Q given as a monotonic decreasing function of an appropriate 
Sg entropy evaluated on /(x). 

We derive the above mentioned entropic upper bound on Q from a special case of the 



Hardy-Littlewood-Sobolev inequality 30| which reads as follows 



Theorem 1. If qi = ^^^^ with < X < D and f e L'^' (M^) then 

f [ f{^)f{y)\^-y\-'d^dy <C{D,X)\\f\\l (16) 



with the constant 
The case of equality is reached if and only if 

2D-A 

/(x)=A(7^ + |x-ap)-^ (18) 
for some A G C, 7 G M\ {0} and a G M^. 

We first remark that this inequality is invariant if x and y are shifted by a so that a can 
be assumed equal to without loss of generality. Function / (x) can then be written as 

/(x) = V(l + ^) ' • (19) 



Since moreover inequality f[T6|) is homogeneous in /, we may also assume that 

[ /(x)rfx=l. (20) 

The probability density f llSp is clearly a g— Gaussian distribution. We remark moreover that 
since < \ < D, 

D 2D-\ 1 

— < = < D 21) 

2 2 l-g2 ^ ^ 



so that 



2 2D-X-2 ^ 1 

1 <q2 = ^ < 1 (22 

D ^ 2D-X D ^ ' 
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The term \\f\\q in the upper-bound of (fT6|) can be recast as, 



CWfWi = C r (x) rfxj = C (1 + (1 - q^)S,,[m , (23) 

which is a monotonically decreasing function of the entropic functional Sg-^ evaluated on the 
density /. 

The integral appearing in the left hand side of equation (|T6l) can be interpreted as pro- 
portional to the binding energy f2 of a mass density distribution p(x) = Mf(x.) (with total 
mass M) where the interaction of mass elements p(xi)(ixi and p(x2)(ix2 separated by a 
distance r = |x2 — Xi| is described by a potential function behaving as 

- -jp(xi)p(x2)cixi(ix2, (24) 

where a is a constant characterizing the "strength" of the above mentioned interaction. The 
alluded binding energy is then 

Q = aM^ [ [ / (x) / (y) |x - y|-^rfxrfy. (25) 

Now, combining equations (fT6|l . (123!) and (l25ll . we get the desired bound on the binding 
energy, 

>(x) 



n < CaM' 1 + (1 - gi) S, 



<?i 



, - (26) 

In the case of the gravitational interaction in D-dimensional space we have X = D — 2. 
Consequently, the entropic parameter qi corresponding to the entropic bound (l26ll is 



In order to get some further physical insight on the particular structure of the entropic 
bound B, it is instructive to consider it from the point of view of its scaling behaviour. It 
is clear that the binding energy Q obeys some simple scaling laws. It seems reasonable that 
a physically meaningful entropic upper bound for Q should satisfy the same scaling laws as 
Q. This requirement is particularly sensible if the bound is exhaustive, in the sense that 
there are density distributions that do achieve the bound. Now, a natural candidate for an 
entropy-related functional of p(x) satisfying simple scaling laws is provided by a power-law 
functional proportional to 
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(28) 

The question is now to find the appropriate values of the exponents q and n. Let's consider 
first the scahng transformation p(x) — )■ p(x) = Ipix), with Z > a scahng parameter (note 
that the mass density p(x) is not in general normalized to 1). The quantities Q and / 
transform respectively as f2 — )■ PQ and / — ?■ /'"^J, which suggests 



p(x)^ dx 



nq = 2. (29) 

On the other hand, consider the transformation p(x) — )■ p(x) = p{l^. In this case Vt and / 
behave as — t- l^^^^Q and / — )■ l^'^^I, suggesting, 

D + 2 = kD. (30) 
Combining f l29|) and fl30|l we finally obtain 

o = , and K = (31) 

which are precisely the values of the exponents appearing in the bound (1261) . This argument, 
of course, does not constitute a formal proof of the bound (the proof, as we already explained, 
is based on the Hardy-Littlewood-Sobolev inequality). However, the above scaling argument 
suggests that it is unlikely that an entropic bound for the gravitational binding energy can 
be found which is not based on the functional (1281) with the exponents ( 13T|) . In other word, 
the entropic bound fl2B]) seems to be unique. 



IV. g-GAUSSIAN MAXIMIZERS 

Let us consider now the particular case of highest relevance for the study of self gravitating 
systems, which is the one corresponding to three dimensional space {D = 3) and A = 1, in 
which case the integral appearing in equation fl25|) ) has the standard form of the gravitating 
binding energy of a mass distribution in three dimensional space. In this case we have 
qi = 6/5 and the entropic upper bound for the gravitational binding energy of a system of 
total mass M described by a density distribution p(x) is. 
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n < i?[p(x)] = C(3,l)GM^(l - i^e/s 



M 



(32) 



Let us now consider the density distributions achieving the entropic bound fl32|) . Com- 
bining equations (fT5|) . (fT9|) . and ( 120|) we can write the maximizer mass distributions as 



p(r) 



3M 

47j-^3 



7^ 



(33) 



where r = |x|. Now, equation ([33 



is well known in astrophysics 



29 



describes the celebrated Plummer density profile, which 



31 



33| . The parameter 7 appearing in f l33p is the "Plum- 



mer radius". It is a scale parameter determining the size of the inner "core" of the mass 
distribution. 

Plummer 's density profile is of considerable theoretical interest. It corresponds to a 
polytropic sphere distribution in position-velocity space with polytropic index n = 5. This 
distribution has the remarkable property of constituting the limit case of a polytropic (and 
a maximum-Sg) distribution exhibiting a finite total mass and therefore being physically 
plausible. As already mentioned, polytropic spheres with n < 5 have finite total mass, while 
those with n > 5 have infinite mass. Another notable feature of the n = 5 polytropic sphere 
is that it is a maximum-Sg distribution in position-velocity space that leads to a density 
distribution (Plummer's distribution) which is itself a g-Gaussian (that is, a maximum-S'^ 
distribution). 



V. CONCLUSIONS 



In the present work we obtained an entropic upper bound B[p] on the gravitational 
binding energy Q associated with a density mass distribution p(x). This upper bound is a 
simple function of an appropriate power-law non-additive Tsallis' entropy Sq evaluated on p. 
In the case of the standard, Newtonain gravitational interaction in three dimensional space 
the alluded entropy is characterized by the entropic parameter q = ^■ 

The density distributions achieving the entropic bound (that is, satisfying Q = B) exhibit 
the form of an isotropic g-Gaussian distribution with <? = §• These maximizer distributions 
coincide with the Plummer density profile, well known in astrophysics. 

It is well known that many of the important facts concerning gravitation are related to 
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simple scaling laws, Kepler's celebrated third law of planetary motion being historically the 



first example 3J]. Here we have provided a heuristic scaling argument suggesting that the 
entropic bound B[p\ for the gravitational binding energy is unique, with the specific value 
of the Tsallis parameter qi appearing in the bound determined by the space dimension D. 

The present findings constitute intriguing new evidence of the relationship that exists 
between the physics of self gravitating systems, on the one hand, and the statistical 
formalism associated with non-additive, power law entropic functionals, on the other hand. 
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